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Abstract 

The statistics of wrinkling flame front is investigated by the quantum field theory 
methods. We dwell on the WKB approximation in the functional integral which is analo- 
gous to the Wyld functional integral in turbulence. The main contribution to statistics is 
due to a coupled field-force configuration. This configuration is related to a kink between 
metastable exact pole solutions of the Syvashinsky equation. These kinks are responsible 
for both the formation of new cusps and the rapid power-law acceleration of the mean 
fiame front. The problem of asymptotic stability of the solutions is discussed. 



1 Introduction 

It had been shown in that under a weakly nonlinear approximation, the dynamics of a 
wrinkled planar flame front is governed by a nonlinear partial differential equation (PDE) 

Here, $ is the interface of a distorted planar flame, Ub is the speed of the planar flame relative 
to the burning gas. Dm is the Markstein diffusivity and 7 is the thermal expansion coefficient, 

7 = (2) 

Pf 

where p/ is the density of the fresh mixture and pb is the density of the burned gas, pf > pb. 
The equation is asymptotically exact in the limit of small 7 ^ 1. A{. . .} represents a linear 
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singular nonlocal operator defined conveniently in terms of the spatial Fourier transform by: 

/ + 00 
dk ^{x,t)e^^'^'' . (3) 

A is responsible for the Darrieus-Landau instability 0, 0]. 

Direct numerical simulations for (JlD performed in ^ show that even when the initial condi- 
tions is chosen to be smooth, the cusps develop on the flame interface as time increases. When 
the integration domain is wide enough, the secondary randomlike subwrinkles arise on the in- 
terface. This wrinkling process is accompanied by the flame speed enhancement undergoing an 
acceleration in time 0], namely the mean radius seems to increase with time according to a 

^3/2 

power-law. 

Numerous analytical investigations devoted to (|l|) display that in the limit of long times 
the local flame dynamics is driven by the large-scale geometry f^-ll^. Exact solutions of (|1|) 
can be obtained in principle by using the pole decomposition technique 0)[|ni]- [|T^]. For such 
pole solutions, (|ip formally reduces to a finite set of ordinary differential equations (ODE) 
which describe the motion of the poles in the complex plane. These poles are interpreted as to 
be related to the cusps observed in physical space. However, numerical and analytical results 
demonstrate convincingly that the solutions of the ODEs do not resemble those obtained from 
the direct numerical integration of (|T]), |]^. In particular, the number of wrinkles obtained from 
the ODEs is independent of time (see [|ll]) and the corresponding (mean) expansion of the 
front is much slower than the t^^"^ power law. 

In [^-[0, it was argued that the inconsistencies with the pole decomposition method lie in 
the stability of the exact pole solutions. The initial value problem of the linearised PDF about 
a pole solution has been solved numerically, as a result they concluded that pole solutions are 
unstable for large 7. Consequently, they are not observed in experiments. 

It was conjectured in that nonlinearity alone is not enough to meet the experimental 
observations and that the results of the spectral numerical integrations is due to computational 
noise. In [^] a model had been developed, where pseudorandom forcing is included. It is shown 
that many broad-banded exciting fields indeed lead to the rapid spawning of wrinkles. 

Analyzing the Fourier spectrum of the solution of PDE obtained numerically in [0, one 
observes that it undergoes the local kinks at very short time. Occurring in the amplitude of 
Fourier modes and its harmonics, these kinks, are followed by an oscillation of this amplitude 
which corresponds, in physical space, to the birth of new pairs of wrinkles. 



The linear stability of the pole decomposition solutions was discussed in |]14[ in details. 
The exact analytical expressions for the eigenvalues and eigenf unctions have been constructed. 
Based on these expressions, in |T^, they demonstrate that for any value of the parameter 7 there 
exists the only asymptotically stable solution with the largest possible (for this particular value 
of 7) number of poles, A^^. As the parameter 7 increases, the equilibrium states of the PDE 
undergo a cascade of bifurcations. In this way the new solution with A^ poles gains stability 
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while the former one with — 1 poles becomes unstable. However, the nonlinear stability and 
dynamics of cusps still remain an important open question within such an approach, 



In the present paper we continue the investigation of the kink type behavior of the expanded 
flame front which has been started in 0-0. In contrast to the study in [|14|, we keep the value 
of 7 fixed and small and use a pseudorandom force as an origin of the spawning of wrinkles. 

We demonstrate that the main contribution to the statistics of wrinkling flame front is given 
by a coupled field-force configurations - the instantons. These configurations are related directly 
to a very short time (practically instant) kinks between metastable groundstates incident to 
different numbers of poles. We should stress that we address the stochastic problem based on 
the Syvashinsky equation (|l|), but not the exact solutions of this equation. 

The paper is organized as follows. In Sec. II, we formulate the stochastic problem for the 
equation governing the slope function dynamics of advancing flame fronts. These fronts usually 
either form fractal objects with contorted and ramified appearance or they wrinkle producing 
self-affine fractals characterized with some critical exponent, ||TB[. In Sec. Ill, we analyse the 
problem from the point of view of critical phenomena theory. We demonstrate that the t^^^ 
power law observed experimentally would be a direct consequence of the general Kolmogorov's 
scaling with the critical dimensions of time Af = —2/3 and velocity = —1/3 (the dimension 
of X is taken by definition as A^. = — 1) which are well known in the fully developed turbulence 
theory (see for a review). 

We must say that in actual problem the renormalization group technique (RG) (which has 
proved itself so well in the fully developed turbulence theory, [l^, [|l^) is ultimately ineffective 



since, obviously, the regime of critical scaling is not attained. One even can hardly use the 
concept of critical dimensions for the actual quantities. 

The examples of successful application of the saddle point calculations to the Burger's 
equation and to the description of intermittency phenomenon in turbulence have been 
given recently. These papers have inspired us to employ this technique in the problem of 
wrinkling flame front expansion. The infinite set of instanton like solutions which we have 
found is dramatically dissimilar to those computed in |^ and p7| . 

In Sec. IV, we construct the statistical theory of wrinkles based on the action functional 
relevant to the actual stochastic problem. 

The minimization of action discussed in Sec. V requires that the field and force to be 
coupled in some particular configurations. We also illustrate the instanton mechanism of poles 
generation for the particular two poles initial configuration. The process keeps repeating itself 
as time increases. We then conclude in the last section. 
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2 The stochastic problem for the Syvashinsky equation 

The stochastic problem for the equation governing the slope function dynamics of the flame 
front, u{x, t) = dx:^{x, t), with a Gaussian distributed pseudorandom force included in the r.h.s. 
reads as follows 

dtu + Ubud.u = DmSIu + ]pUb^{u} + /. (4) 
The pair correlation function for / is taken in the form 

(/(x, t)/(x', t')) = Df{x- x')6it - 0, (5) 

in which the function Df{x — x') is supposed to be an even smooth "bell" -shaped function of 
X. To be specific, we take it in the form 

n r ^ Dq rri 

Df[x) = — (6) 

■K x^ + 

decaying at the rate m and turning into Do6{x) as m — 0, where -Do is a constant. 

The equation @ is similar to the Burger's equation except the singular term, oc ■yUbA{u}. 



The homogeneous unforced equation (^) considered in [11| in details. In particular, it was 
shown that it possesses a pole decomposition, i.e., it allows a countable number of uniform 
solutions, 

in which z'-s are poles in the complex plane (coming in complex conjugate pairs) moving ac- 
cording to the laws of motion of poles 

Zi = -2u — ^ i'^Ui, sign [^{zi)] , (8) 

where denotes the imaginary part of a pole. One can derive easily the corresponding steady 
{zi = 0), solution of the Sivashinsky equation for the simplest configurations concerning the 
minimal number of poles. For example, for two poles the only steady solution is given by 



ADmX 



= -^fT^^ (9) 

and there are two possible four-pole steady configurations, 

^(4). X _ ^ 4Dm(±2x3 + 2lV2iDl, + 9V2iDmx^) 
" -x^±hAy/2iDljX±Q^iDMX^ + SlDlj' 

In ||5|-|0], by numerically solving the initial value problem that results from linearizing the PDE 
about a pole solution, it was concluded that such pole solutions with a fixed number of poles 
are asymptotically unstable, while the actual solutions which can be observed experimentally 
undergo a sequence of kinks between different metastable groundstates of the type (0). 
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To construct the solutions with spawning wrinkles, we exploit the exact correspondence 
between an arbitrary stochastic dynamical problem with the Gaussian distributed random 



force and a quantum field theory, [^. A short and elegant proof had been given in p2|. The 
stochastic problem (^) and (|^) is completely equivalent to the field theory of two fields with an 
action functional 

S[u, ^] ~ 2 J '^^ ^-^ '^y ^{^^ t)Df{x — y, t)w{y, t) (11) 

Here, w{x,t) is the auxiliary field which comes into play instead of the random force /. w{x,t) 
determines the response functions of the system, for instance, the linear response function is 
{uw) . 



3 The power law of the flame front expansion 

The mean squared distance of propagating fiame front, B?{t), can be expressed naturally via 
the linear response function mentioned at the end of the previous section as following, 

R^{t) = J dx x'^ {u{x,t)w{x,0)) . (12) 

The requirement that each term of the action functional be dimensionless (with respect to x 
and t separately) leads to the power counting relation for the product, A[mw] = d, where d is 
the space dimensionality. Therefore, the power of the linear response function is 

A[{uw)] = d- d = 0. (13) 

We note that (|TB|) is still valid whether a critical regime is attained or not. 

Following a tradition, we accept the natural normalization condition that A[x] = —1, then 

A[i?2] = -2. (14) 

From the other hand, the observation data show that 

oc t^ (15) 

(i.e., R oc t^/^). Comparing ( [T^ ) and (P^D, one concludes that, in the theory (|TT|), the time 
variable t possesses the effective dimension, 

A[t] = -2/3, (16) 

which is equal to the Kolmogorov's critical dimension of time in fully developed turbulence. 

If the result ([16D were a true critical dimension, then one could claim that there is a criti- 
cal regime in the wrinkling flame front propagation problem, i.e., for any correlation function 
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there is a definite stable large scale long time asymptotics. Nevertheless, we stress the dra- 
matic difference between the stochastic theory of turbulence (the Navier-Stocks equation with 



a random forcing included, [0) and the actual problem. From the point of view of the critical 
phenomena theory, the problem - is formulated erroneously. 

Namely, interesting in the long time large scale asymptotics behavior of correlation functions, 
one has to omit the term oc Dm^'^ (in the momentum Fourier space) from the action ([TT| ) in 
benefit to oc '~fUb\k\. However, in this case, there are infinitely many Green's functions which 
have singularities with respect to a general dilatation of variables, i.e., such a theory cannot be 
renormalized. 

One can investigate a theory in which the both terms are included simultaneously. Up to 
our knowledge, a model where the concurrence between two terms (in the momentum Fourier 
space) oc A;^ and oc |A;p~^° (0 < a < 1/2) was considered firstly in in the framework of 
the renormalization group approach. It is shown that up to the value ac < 1/2 a regular 
expansion in a and e (the deviation of the space dimensionality from its logarithmical value) 
can be constructed and then summed over by the standard renormalization group procedure. 
The critical indices of all quantities are still fixed on their kolmogorovian values. 

However, for a > 1/2 and for the particular case of = 1/2 which we are interested in, 
the renormalization group method fails. The matter is in new additional singularities which 
spawning in the correlation functions of the field u as \k\ 0. Such singularities cannot be 
handled by the renormalization group in principle since they do not related to a general scaling 
with respect to dilatation of variables. 

The summation of the leading infrared (|A;| 0) singularities of correlation functions can 
be done by an infrared perturbation theory. If we limit ourselves to the functions Uz{x, t) which 
have poles {z{t)} in the complex plane, then the action of the singular operator A is reduced 
to a first order derivative operator 

A{uz} = i sign C^lzit)]) d^u^ix, t), (17) 



see 



11 1 . Then, in the momentum-frequency representation, the term with A can be taken into 



account ClS db small shift of frequency (oc 7), 

u ^ u - i-fUb\k\. (18) 

The infrared perturbation theory results from the expansion of exp S over nonlinearities. The 
corresponding diagram technique coincide with the diagram technique of Wyld, . The lines 
in the diagrams are associated with the bare propagators, in the Fourier space. 

Gun. = Gl, = + (19) 

and 
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where Df{k) is the momentum representation of (|^). For any correlation function, this diagram 
technique gives an infrared representation which is naturally consistent with the 7 expansion 
and is well defined for small values of the parameter 7. 

We are not going to discuss the application of the infrared perturbation theory to the 
actual problem in detail. Here, we conclude that for any pseudo-differential operator oc 
< a < 1/2, the model of the type ([TTl) has a critical regime with the critical indices fixed at 
their Kolmogorov's values (see @). However, for a > 1/2 the stability of asymptotics still an 
important open question if 7 is large. 



4 Statistics of planar flame front wrinkling 



We are going to discuss the saddle point configurations of ([TTD which can provide us by a detailed 
description of mechanism of wrinkles generation on the propagating flame front surface. For 
the future purposes, it would be convenient to perform consequently the rescaling of fields in 

(0), 

u 

—, w ^ UbW, (21) 

such that the parameter Ub is removed from the nonlinear term wuUx, and then another rescal- 
ing, 

u^-, Ub^^. (22) 

7 7 

As a result of such a simple transformation, we observe that the parameter of thermal expansion 
7 which we assume to be small, plays the formal role of h in quantum field theory: 

S -S. (23) 

7 

The correlation functions of the basic field u are then given by the functional integral 

Gn{,Xi,ti] X2, t2]... Xn, t^) = ( Vu Vw m(Xi, ti)u{x2, ^2) • • • u{Xn, t^) exp(--S'), (24) 

J 7 

and can be derived naturally by means of a generating functional which has been introduced 
firstly in |2^ and then employed in ||2^, |2^ 



Z{X) = (exp J dt dx -^""^^ = J 6xp ^—{—5' + i J dt dx Am}^ . 



(25) 



The coefficients of the expansion of ^ in A are the correlation functions (U). 

There are no general methods to compute such a functional integral exactly. The straight- 
forward perturbative approach is to expand the exponential in the functional integral (^) in 
powers of the nonlinear term wuUx- However, since we are interested in nonperturbative effects, 
it seems more natural to search for some saddle-point configurations which minimize the action 
functional ([TI|), thus dominating the functional integral in a way similar to the saddle point 
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approximation in ordinary integrals. Such solutions are called instantons, and they determine 
the asymptotics of (p5D at small 7 <^ 1 which corresponds to WKB approximation in quantum 
field theory {h <ti 1). 

Another quantity which can be expressed via the generating functional (^) is the probability 
distribution function V{u) for the field u, 

V{u) = y Z{\) exp (^-i j dt dx \v}j . (26) 

The behavior of V{u) for large u is also dominated by some saddle-point configurations of the 
integrand. However, these configurations are not the same for both (^) and (^Gl) . 



5 Kink solutions of the forced Syvashinsky equation 

In what follows we shall look for saddle point configurations driven by the random force in 
terms of functions which have poles in the complex plane, Uz{x,t). To be specific, we observe 
generation of the four poles configuration from the two poles as a result of a kink. In contrast 
with 1^ and we need not introduce here a large artificial parameter to fix the saddle points 



dominating the functional integrals (^) and (|26D since we have the inverse thermal expansion 
coefficient I/7 which is naturally large. 

Suggesting that the field u can be continued analytically on the complex plane except for 
the poles, we shall study the correlation function of the form 



G{z) = ( exp 



u{z) - u{z*) 
7 



(27) 



of two distinct points of the complex plane symmetrical with respect to the real axis. We 
suppose also that at the initial moment of time t = the field u can be depicted as a config- 
uration of two complex conjugated poles, z and z* . The function (^) possesses a generating 
property: Tailoring (pTj) in powers of I/7, one obtains the "structure functions" for the field u. 
The functional Fourier transform ( p6D of ( p7[ ) gives us the two point probability distribution. 
The structure function generated by (E^ is related to the same point x = ^{z) on the real axis. 



Taking an average in (|27|) with respect to a functional measure, we perform an integration 
over all possible configurations u{x, t) with the asymptote prescribed by initial two poles and all 
possible final multipole configurations. The basic symmetry of the action (|1^) is the Galilean 
invariance which reveals itself in the real transformation 

Ua{x,t)^u{x + Xa{t),t)-a{t), (28) 

where a{t) is an arbitrary function of t decreasing rapidly as \t\ oo, and Xa{t) = Jq dt' a(t'). 
The transformation (|28|) defines an orbit in the functional space of u along which the result of 
functional averaging does not change. It follows that the integral itself is proportional to the 
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volume of this orbit. This volume should be factorized before one can perform the saddle-point 
calculation (see [^]). It is appropriate to choose for the latter the "plane" transversal to the 
real axis ^{z) and then cancel out the real components of u which are related to each other 
via ([28|). In (^), the real contribution to u is subtracted out, so that it is very suitable for 
instanton calculations. 



The asymptotics of small 7 in (|27|) is dominated by the saddle point configurations of the 
functional 

W[u,w,z] = ""(^^ - _ s[u,w], (29) 
7 



which should satisfy the following equations obtained by varying of (^9]) with respect to u and 



ut + uux - Dmu^x - \UhJ^{u} = -\— [ dx' Df{x- x', t)w{y, t), (30) 

/ 2 •y J 

1 i 

wt + uw^ + Dmw^x + -UbA{w} = — 6(t) {6{x - z) - 6{x - z*)} . (31) 

2 7 



These equations for the saddle point configurations are similar to those derived in [26| except 
the last singular term in the l.h.s. They follow from the Syvashinsky equation for the slope 
function (Q), however they contain information on a special force configuration necessary to 
produce instantons also. 

Indeed, the particular solutions of (^) and ( pID are dependent substantially from the initial 
data for u and w. Minimization of the action requires m — > 0, at t ^ —00 and w —>■ 0, at t ^ 00. 



Obviously, any solution of the equation (0) which is nonsingular as t ^ +00 should be equal to 
zero at t > (since the field w feels a negative diffusivity) . Following an analogy with p6| , ^], 
one can say that the field w propagates backwards in time starting from its initial value 

w{t = -0) = -— {6{x - z{0)) - 6{x - z*{0))} (32) 
70 

while it is zero at all later moments of time. Therefore, the system (|3D|) and (^Tj) as well as the 
integrals in (|TT|) can be treated for t < only. 

While propagating backward in time, (^) is a subject to a drift of the initial conditions as 
governed by the velocity in the equation (|3TD, the smearing of the initial 5— function distribu- 
tions in ( p2D due to diffusivity, and finally, an advection in the complex plane, in the imaginary 
direction towards the real axis. In the limit of no diffusivity, one can neglect the smearing in 
the equation (|3lD . A simplified equation, which we arrive at when we drop the diffusivity term 
is just moving the 5— singular right hand side of (0) around. Therefore, the solution of (pT|) 
can be expressed naturally in the form 

""^^^ ^ ~W) ^^^"^ ' '^^^^ ~ ^^"^ " ^^^^ 

with the boundary conditions 7(0) = 70, z{^) = zq. If one takes the diffusivity term in (^) 
into account (in the case of eventually small Duj^Ub < 1); the solution of (p^) would be 



Instanton solutions in the problem 



10 



expressed in terms of even functions decaying as |x| ^ oo, 

^(t) = '—I yM \ (34) 

where the parameter 

^« ^ '^^> 

is a size of an instanton, and z{t) is its position changing with time (we have borrowed the 
terminology from the quantum field theory). As \y{t)\ 0, the instanton shrinks to a point, 
and the solution ( PH ) is reduced to ([331). 

To proceed with the equation (pOf), we rewrite it due to ( p!7D in the form 

l^Ub 

7 

The one-dimensional equation (|36|) can be linearised by the Cole-Hopf transformation, 



i 1 iUh f 

ut + uu^ - -Da/Mxx - -Ubsign '^[z{t)]d^Uz{x, t) = — / dx' Df{x - x', t)w{y, t). (36) 



so that we arrive at the equation 



tpt - Dm^^xx - -UbSign Q[z{t)]d^tlj{x,t) = — / dy D f{x - y,t)w{y,t). (38) 
2 2 ^0 J 

Since w{x,t) distinguishes from zero only for t < and Df oc 6(t), the only nontrivial contri- 
bution into the r.h.s. of ( pS] ) is given by the moment of time t = 0. The general solution of the 
equation ( |55| ) in the Fourier space reads as follows 

i^^nstik,t) = l-—Df{k,0)w{k,0)6{t), (39) 
^ 7o 

plus a transient process decaying rapidly as time growing, oc exp [— t(Djvf^^ + f^fol^l)] • 



Now we can use (y) and (|3^) to write down the r.h.s. of (|39|) explicitly, 

Df{k,0)w{k,0) = A7r^iDoexp[-{y{0) + m)\k\ - ik^zo]]smkQ[zo]. (40) 



Performing an inverse Fourier transform of (p9l), one obtains 



inst 



27r 



{Dm + l^U.mf + 72f/2(x - ^o)^ (^M + 7of/6m)2 + ^2f;2(^ _ ^*)2 

" (41) 

Finally, we arrive (||29|) at the following four poles configuration for the instanton Uinsu 



ADmx^\z] 

Uinst = S(t). 42 

{{DM/-foUb + m)2 + x^- Q[z]^f + Ax^Q[z]^ ^ ^ ^ ^ 

For the last step of instanton computation we have to define the functions 7(t) and ip{t) = 
'^[z{t)] in the equation (Rl) (remember that ^[z] is fixed) as t < 0. We can do it by a direct 
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substitution of ( P^D (in case of eventually small diffusivity) into the equation (^). Here we 
note that since the instanton solution (|4^ ) exists as t > (if one takes the transient process 
into account), we have to use the initial two pole configuration instead of u{x,t) in (|3TD. As a 
result, we obtain the system of simplified equations 

4:xDm 



at) 



The formal solution of ( ^3D is given by 

7(t) = 7oexp - 



a;2 + if{ty 
Ax Dm 

x2 + ifity 



AxDm 



lit), 

:dt'] , 



it' < 0) 



It x^ + ip{t'y J 

and can be computed, in principle, if one knows '^{t\ The equation 



is equivalent to 



X 



D 



M 



lnif{t) + 



which leads to 



2xDm 
tDM 



X 



^p{t) = v?oexp 

W{x) is the Lambert function which meets the equation 

W{x) exp W{x 



-2tDM/x' 



X^ 



X. 



(43) 
(44) 

(45) 

(46) 
(47) 

(48) 



The latter equation has an infinite number of solutions for each (non-zero) value of x. W has 
an infinite number of branches numbered by an integer number n G [— oo . . . oo]. Exactly one 
of these branches is analytic at (the principal branch, n = 0). The other branches all have a 
branch point at 0. The principal branch is real- valued for x in the range — exp(— 1) . . . oo, while 
the image of — oo . . . — exp(— 1) under W{x) is the curve —ycot{y) +yi, for y G [0 . . . vr]. For all 
the branches other than the principal branch, the branch cut dividing them is the negative real 
axis. The image of the negative real axis under the branch W{n, x) is the curve —y cot{y) + yi, 
for y G [2kTr . . . {2k + l)rr if k > and y G [{2k + l)7r . . . {2k + 2)n] if k < -1. These curves, 
therefore, bound the ranges of the branches of W, and in each case, the upper boundary of the 
region is included in the range of the corresponding branch. 

Each particular orbit of (|47|) provides a distinct solution of (|45|) and (p3|). However, each 
configuration w{x, t), t < which enjoys (^) is related to the same configuration u{x, t),t>0. 
The value of (^) for the instanton is obviously finite, however, one can hardly compute it for 
each branch of w{z,t). 

Let us consider the principle branch of the function ip{t) just to illustrate the idea of com- 
putation. One can check that the leading contribution to if is accumulated around x = 0. The 
asymptotic behavior of W at complex infinity and at is given by 

log(log(x))('"+^) 
log(x)(™+"+i) ' 



W{x) ~ log(x) - log(log(x)) + ^ C( 



(49) 
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where log(x) denotes the principal branch of the logarithm, and the coefficients C{m,n) are 
constants. 

Restricting to the ffist term of the asymptote (|49|), one obtains 

ip{t) - -ip{0)DMt, {t < 0), (50) 

then we use (0) and (|i5| ) to compute 7(t) as t < 0, 

7~7oexp — tan , t < 0. (51) 

Now it is a matter of a simple computation to find the action on the principal branch of the 
instanton. We collect everything together and substitute (0), (|50|), (^2]), and (^) back to (0) 
to obtain 

'S'mst = ^exp[ tan^Vo], (52) 

m ipo 

while the correlation function we have been studying is 

G oc exp f — ^exp[ tan""*^ o^ol | . (53) 

\m ipo J 

This simplified formula is obviously not an exact answer. It is just a leading asymptotic of 
G'(53[2;o]) if 1/70 is a large number, ipQ = ^'[^^o] is eventually small, however, not too small (since 
we have not taken the smearing due to diffusivity into account). Nevertheless, for some interval 
of scales '^[zq], rather small than large (see Fig. 1), the asymptotics prescribed by (^3|) (a thick 
line) is very close to a model curve which corresponds to the Kolmogorov's critical exponent 
for velocity, A^[f] = — 1/3. . 

Due to a specific property of the action (|11]), the contribution from fluctuations of 6w and 
6u up to second order will be zero. To find higher order corrections to (|55D, one has to consider 
at least the third order terms. We investigate this problem in future publications. 



6 Discussion and Conclusions 

The study performed in this paper conffims that the stochastic model for the outward propagat- 
ing flame in the regime of well developed hydrodynamic instability leads to the self-fractalization 
of the flame front. The mechanism consists of successive instabilities through which the inter- 
face becomes more and more wrinkled as time increases. The main contribution to the statistics 
of wrinkles is due to a coupled field-force configurations which are found to be responsible for 
the birth and growth of wrinkles. This behavior is different from exact pole solutions for which 
the number of poles is constant. The acceleration of the mean front radius is clearly due to 
successive births of poles. 

We have demonstrated that the t^^'^ power law observed experimentally for the mean radius 
of the advancing flame front would be a direct consequence of the Kolmogorov's scaling with the 
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critical dimensions of time Af = —2/3 which is well known in the fully developed turbulence 
theory. Provided the critical regime in the model of wrinkling flame front exists, then it 
means that each correlation function in the model has a definite stable long time large scale 
asymptotics, which does not depend from the particular sequence of kinks. If one replaces the 
nonlocal operator A with some pseudo-differential operator oc < a < 1/2, then, as 



it was shown in |£4|, the model of the type (|TT|) has a critical regime with the critical indices 
fixed at their Kolmogorov's values. However, for a > 1/2 the stability of asymptotics still an 
important open question if 7 is large. 

We have used the saddle point calculations assuming the inverse thermal expansion co- 
efficient 1/7 as a large parameter. As a result, we construct an infinite family of instanton 
solutions numbered by n G [—00 ... 00]. Each instanton determined by one of the branches of 
the Lambert function W{x) has a unique behavior as t < , however, all instantons are undis- 
tinguishable as t > 0. The asymptotic behavior of these solutions is close to that of prescribed 
by the Kolmogorov's scaling. 

The crucial problem of the developed technique is of contribution to the action from the 
fluctuations against the instanton background. The general analysis of the set of instantons 
which we have found will be published in a forthcoming paper. 
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